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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering 
Development  Center  (AEDC),  Air  Force  Systems  Command  (AF5C),  under 
Program  Element  65807F.  The  results  were  obtained  by  ARO,  Inc.,  AEDC 
Division  (a  Sverdrup  Corporation  Company) , operating  contractor  for 
the  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee.  The  work  was 
done  under  ARO  Project  No.  R33A-03A.  The  author  of  this  report  is 
J'.  C.  Chlen,  ARO,  Inc.  Elton  R.  Thompson  is  the  Air  Force  project 
manager.  The  manuscript  (ARO  Control  Ho.  ARO-ETF-TR-77-44)  was 
submitted  for  publication  on  June  17,  1977. 
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1.0  INTRODUCTION 

Conical  diffusers  are  commonly  used  In  wind  tunnels,  test  facili- 
ties, propulsion  systems,  etc.  As  the  flow  moves  through  the  diffuser, 
the  velocity  decreases  and  the  static  pressure  increases.  The  amount  of 
the  pressure  rise,  l.e.,  the  diffuser  performance,  is  strongly  influenced 
by  the  diffuser  inlet  flow  conditions  and  the  divergence  angle.  In 
addition,  the  inlet  boundary- layer  thickness,  the  velocity  profile,  the 
inlet  Mach  number,  Reynolds  number,  and  the  turbulence  intensity  all 
affect  the  diffuser  performance  in  such  a complicated  way  that  diffuser 
design  has  relied  almost  solely  on  empirical  data. 

Using  recent  advances  in  computational  fluid  dynamics  and  turbu- 
lence modeling,  a computational  method  has  been  developed  for  incom- 
pressible diffuser  flows  (Ref.  1).  The  method  allows  one  to  obtain 
numerical  solutions  of  the  Navier-Stokes  equations  in  finite-difference 
form.  A two-equation  k-e  model  is  employed  for  the  turbulent  flow  along 
with  a sublayer  coordinate  transformation  so  that  the  flow  field  for  the 
entire  diffuser  can  be  computed  for  both  separated  and  nonseparated  flow 
cases . 

The  purpose  of  the  present  study  Is  to  further  verify  the  numerical 
method  developed  in  Ref.  1 by  applying  it  to  a conical  diffuser  configu- 
ration with  various  inlet  boundary-layer  conditions.  The  development  of 
a two-equation  turbulence  model  suitable  for  axisymraetric  flows,  a 
multiple  grid  Gauss-Seidel  iteration  procedure  for  faster  convergence, 
and  a series  of  numerical  computations  are  presented  along  with  compari- 
sons with  available  experimental  data. 

2.0  GOVERNING  EQUATIONS 

The  basic  governing  equations  for  the  steady-state  incompressible 
turbulent  flow  can  be  written  as: 
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Continuity  Equation: 

5 (r5„)  + v)  = 0 

OX  OT 


CD 


X-Momentum  Equation: 


§?)]} 


(2) 


Y-Momentum  Equation: 


dx  <3r  dr  \p)  \ dr  ^ dr  -r  5x  1_  \dx  dr/j 


+ 2 


^(fr  - f)} -{(!)'■  5^ 


(3) 


The  turbulent  quantities  such  as  u’v'  can  be  written  in  terms  of 
the  mean  flow  quantities  through  the  eddy  viscosity  concept  as 

'*  ^ ~ I*' 


v'2  = -2.,  |k 

‘ dr  3 

where  is  the  eddy  viscosity  and  k is  the  turbulent  kinetic  energy 
(TKE) . 


(4) 


After  substituting  the  eddy  viscosity  expressions  Into  the  momentum 
equations  (Eqs.  2 and  3),  one  can  obtain  the  vorticity  equation  through 
cross  differentiation.  If  one  also  Introduces  the  stream  function 
concept  into  the  continuity  equation,  the  governing  equations  can  be 
replaced  by  the  vorticity- stream  function  formulation  as: 
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Vorticity  Equation: 


(5) 


where 


Stream  Function  Equation: 


where 


(6) 


A second  order  pressure  equation  can  also  be  derived,  as  shown  below. 


idx"  2 ) V 2 ) 


(7) 
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3.0  A TWO-EQUATION,  LOW-REYNOLDS-NUMBER  k-e 
TURBULENCE  MODEL 


The  eddy  viscosity,  is  calculated  through  the  Prandtl-Kolmogorov 
formula,  l.e., 

1.2 


= V • - 


where  the  turbulent  kinetic  energy,  k,  and  its  dissipation,  e,  are 
obtained  from  the  partial  differntial  euqatlons  presented  below. 

Turbulent  Kinetic  Energy  Equation  (k) ; 

^ i - . 1 )L-^l3k  r_^:^  - g + dk  I 

* ^ (^1  "0/  *]  * (if  * 1;)'} 

■ . ^ r y(''A 


Turbulent  Kinetic  Energy  Dissipation  Equation 


" " (if)' 
* (t1  *]  * (If  ’ If)' } 


- C 


2 iv  + v/a^)  k 


T-  = 0 


The  coefficients  used  in  the  present  analysis  are 


c = t + Ac 
fi  fio  (I 


A 


(t°  3 (a  + A/b) 


(8) 


(9) 


(10) 


(11) 
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where 


Cj  = 1.36  . t'2  = 1.92  [l  -0.3  exp.  (-R^)] 

= 1.0 

A = p - H = 

a = 1,100  , b = 0.27 


(12) 


The  coefficient  c^^  is  related  to  the  two-dimensional  (2D)  planar 
flows  while  the  coefficient  Ac^  is  attributed  to  the  axisymmetry  of  the 
flow.  The  coefficient  can  also  be  written  as 


Ac  = ^ 


(13) 


through  the  use  of  the  continuity  equation.  Justification  for  the  form 
of  Eqs.  (11)  and  (13)  is  given  In  the  next  section. 


4.0  DEVELOPMENT  OF  THE  EDDY  VISCOSITY  COEFFICIENT,  c», 
FOR  AXISYMMETRIC  FLOWS 


The  eddy  viscosity  coefficient,  c^,  developed  in  Ref.  1 for  2D 
flows  is 

= A 

fio  3(H+.\/h)  (14) 

where  A,  a,  and  b are  given  in  Eq,  (12).  The  asymptotic  value  of  c^^ 
for  large  A (i.e.,  away  from  the  wall)  is  0.09.  When  the  model  given  by 
Eq . (14)  Is  used  for  conical  diffuser  flow  calculations,  it  produces  a 
faster  centerline  velocity  decay  than  that  indicated  by  experiment  (Ref, 
1).  A similar  result  was  observed  in  the  case  of  a round  Jet  (axisjon- 
metrlc)  with  a two-equation  k-e  model  developed  by  Launder,  Rodl,  and 
Spalding  (Refs.  2 and  3). 
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In  order  to  study  the  effect  of  on  the  flow-field  solution  for 
a diffuser,  was  lowered  to  a value  of  (0.56  throughout  the  flow 

field,  where  is  given  by  Eq.  (14).  The  resulting  centerline  velocity 
is  shown  in  Fig,  1 , It  can  be  seen  that  the  calculated  results  using 
the  lower  value  of  c^  compare  well  with  the  experimental  data  by  Okwoubl 
and  Azad  (Ref.  4).  However,  the  approach  of  simply  setting  c^  = 0-56 
c^^  is  not  completely  acceptable  for  diffuser  flows  because  it  changes 
not  only  the  flow-field  solution  but  also  the  inlet  condition. 

In  the  study  of  axlsymmetrlc  jet  flows.  Launder,  et  al. , (Ref.  2) 
designed  a formula  for  c^  so  that  the  round  jet  flow  field  can  be  computed 
reasonably  well  with  a two-equation  k-e  model.  The  model,  based  on  the 
Jet  centerline  decay,  is 

= 0.09-0.04f  (15) 

where 


YC 

A-  1 

2Au 

V dx  dx  1/  1 

The  exponent  of  0.2  in  Eq.  (16)  is  derived  from  numerical  optimization 
for  the  jet  flow  computation.  The  main  objections  to  this  model  are; 
first,  the  form  of  c is  Somewhat  arbitrary,  and  second,  c in  Eqs,  (15) 
and  (16)  is  a function  of  x only,  i.e.,  c^  Is  constant  in  the  radial 
direction.  Thus,  a new  model  of  c^  was  derived  for  conical  diffuser 
flows.  Note  that  the  new  model  must  also  be  applicable  to  axisymmetric 
jet  flows  because  a jet-like  flow  can  be  obtained  in  a diffuser  by 
increasing  the  diffuser  divergence  angle.  The  requirements  for  the  new 
model  for  c^  are:  (1)  it  must  be  a simple  function  of  the  flow-field 

variables,  (2)  it  must  be  a function  of  both  radial  (r)  and  axial  (x) 
coordinates,  and  (3)  it  must  recover  the  2D  planar  form  when  the  flow  is 
locally  two  dimensional. 
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4.1  A NEW  Cfi  MODEL  FOR  AXISYMMETRIC  FLOWS 

The  shear  stress  can  be  written  in  terms  of  the  velocity  gradient 
in  the  Prandtl-Kolmogorov  model  as 


' \dy  dx  j 

f Jii  ■*-  ^ 

It  f ' dJ 


(17) 


Let 


c = f - Ac 

H fio  ft 


(18) 


then 


- = (c  +Ac^— + 
p \ftc  e \d,  dJ 


(19) 


where  Ac^  is  the  correction  term  for  the  axisymmetric  flows. 

Since  the  c^^  model  works  well  for  both  planar  and  axisymmetric 
pipe  flow,  it  can  be  expected  that  Ac^  becomes  Important  when  3u/3x  is 
nonzero.  To  retain  the  two-dimensional  applications  of  the  model,  one 
can  relate  Ac^  to  3u/3x  + 3v/3y  which  vanishes  in  planar  cases.  Thus 
for  axisymmetric  flow. 


Ac,  - + S 

Ux  dt) 


(20) 


The  above  relation  can  be  written  in  the  dimensionless  form  by 
Introducing  a local  velocity  scale  and  a local  length  scale: 


Ac  = - c , ^ ^ 

ft  (ll  T 

m 


(21) 


where  c . is  a constant,  J.  is  a turbulent  scale,  and  Au  is  the  charac- 
U I K m 

teristic  velocity.  By  equating  the  length  scale  A„  to  and  taking 

lx 


Au  as  u , Eq,  (21)  becomes 
m c 


Ac  = -c  , ^ . aE 

ft  ftl  r f 


(22) 
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The  coefficient  c^^  could  be  determined  from  known  distributions  of 
^c’  using  Eq.  (21),  However,  since  detailed  experimental  data 

for  the  conical  diffuser  are  not  available,  the  constant  c^^  will  be 
determined  from  round  free- jet  data  and  verified  for  conical  diffuser 
flows.  At  the  centerline  of  a round  jet,  v/r  is  0.042  5) 

and  the  length  scale,  is  1.2  Therefore,  one  can  write  Ac^  as 


0.0+2  u 

Ac  = -c  , ^ 

/X  r. 


1.2  t 


iZi  = 


1/2 


Ml 


. (0.05) 


(23) 


Rodi's  estimate  of  c^  at  the  centerline  of  a round  turbulent  jet  is 
about  0.04  (Ref.  3).  Therefore,  from  Eq.  (23) 


c 


Ml 


0.03 


c,  - 0.0+ 
M« 

0.05 


_0.09  - 0,04  = -i  0 
0.05 


(24) 


Thus,  the  final  expression  for  the  eddy  viscosity  coefficient  c becomes 


c = 
M 


c + Ac 
M"  M 


(25) 


where 


jio 


A 

3 (a  + A/b) 


Ac 

M 


= -S  ■ 


y . k , 'fir 
r ' e ' u 

C 


4.2  THE  Cfi  DISTRIBUTION  FOR  A ROUND  JET 


The  derived  c^  applies  to  the  round  jet  flow  with  the  mean  radial 
velocity  distribution  given  by  Schlichting ' s and  Tollmlen’s  theories. 
The  predicted  c^  distribution  Is  compared  with  Rodi's  estimated  distri- 
bution that  was  obtained  from  experimental  data. 


Assuming  the  length  scale  variation  is  small  across  a round  jet, 
1*2  ^2_/2  that  c^^  is  0.09,  becomes  a function  of 

(v/r)  only.  From  Eq.  (25), 


V ■ -(f)  ■ ^ 


(26) 


(v/  u 

0.09  - 1.2  , 

(r/ri^2 
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The  expression  v/u^  can  be  obtained  from  either  Schlichting's 
(Ref.  6)  or  Tollmlen's  theory  (Ref.  7).  From  Schlichting's  theory. 


where 


1 

4 


J 

Bit 


uilh 


''k 


0.0161 


(27) 


Thus, 


(28) 


Also,  from  Ref.  (6), 


and 


^1/2  = = (5.27  * 0.0161)  x 


K 


\/2  (5.27  ♦ 0.0161) 


X 


4 \'7T 


5.27  . vTs 


— ■ 1 


(29) 


Substituting  Eqs. 

c becomes 

y 


(28)  and  (29)  into  Eq.  (26),  the  final  expression  for 


= 0.09  - 0.0424 


i'-W) 


(30) 


The  distribution  of  Eq,  (30)  is  shown  in  Fig.  2.  For  large  values  of 
n,  i.e,,  at  the  outer  edge  of  the  jet. 


c = 0.09+  M626  . 

IL  ^2  (31 

Therefore,  the  asympototlc  value  at  n ^ “ is  0.09,  This  is  consistent 
with  the  flow  behavior,  because  as  one  moves  away  from  the  centerline 
the  factor  (1/r)  diminishes  and  the  flow  behaves  more  or  less  like  a 
two-dimensional  planar  flow.  The  distribution  obtained  with  Toll- 
mien's  velocity  distribution  is  also  shown  in  Fig.  2.  Rodl's  estimate 
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of  for  a round  jet,  based  on  Wygnanskl,  et  al, , measurements,  is  also 
Included  for  comparison.  Although  the  trend  and  level  of  c^  near  the 
centerline  region  is  roughly  the  same,  the  behavior  of  Rodl’s  is  quite 
different  from  the  present  theory. 

5.0  BOUNDARY  CONDITIONS  AND  COORDINATE  TRANSFORMATION 

5.1  BOUNDARY  CONDITIONS 

Boundary  conditions  for  the  diffuser  flow  can  be  divided  into  four 
regions,  namely,  the  inlet  condition,  the  exit  condition,  the  symmetry 
condition,  and  the  wall  condition. 

5.1.1  Inlet  Boundary  Condition 

Two  types  of  inlet  conditions  are  used  in  the  present  analysis, 
the  fully  developed  pipe  flow  inlet  condition  and  a boundary-layer  inlet 
condition  with  a potential  core  region. 

For  the  fully  developed  pipe  flow  condition,  the  profile  is  ob- 
tained by  solving  the  three  stations  in  the  inlet  section  of  the  dif- 
fuser until  convergence  is  reached.  The  convergent  pipe  flow  solution 
is  then  held  fixed  as  the  inlet  condition  for  the  subsequent  diffuser 
flow- field  computations.  With  this  condition,  the  only  parameter  needed 
to  prescribe  the  profile  is  the  inlet  Reynolds  number;  all  other  initial 
condition  parameters,  i.e,,  velocity  profiles,  TBCE  profile,  etc.,  are 
computed  and  are  self-consistent.  For  a nonfully  developed  profile, 
i.e.,  a boundary  layer  with  a potential  core  region,  prescription  of  the 
inlet  condition  is  more  involved.  Every  flow  parameter  (J2,  u,  v,  'F, 
k,  e,  and  v^)  must  be  given  consistently,  not  only  in  the  core  region 
but  also  throughout  the  sublayer.  In  the  present  analysis,  the  boundary- 
layer  thickness  (y^)  is  first  given  and  the  total  shear  stress  distribu- 
tion is  calculated  from 

~ jl -(-  cos  •0<y<Vg  (Boundary  Layer) 

/-  ' “ (32) 

’^ipt  ~ 0 ’ >'>ye  (Core  Region) 
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The  mixing  length  distribution  is  obtained  as 


im  = 0.09  tnnh  /'Mi 

\0.09  - 


(33) 


for  both  the  boundary  layer  and  core  regions.  Equations  (32)  and  (33) 
are  then  used  to  compute  the  vortlcity  distribution  through  the  modified 
Van  Driest  formula 


where 


n+  = ^ - 

i?y+ 

1 

(34) 


The  velocity  and  the  stream  function  are  then  obtained  from  Eq, 
(34)  through  a numerical  Integration  scheme.  The  free-stream  values  of 
the  turbulent  kinetic  energy  and  the  eddy  viscosity  are  computed  from 


k = 1 (kl 
0 2 


=12 


j 


"I.C 


(I) 


(35) 


where  kl  and  kZ  are  free-stream  proportionality  constants.  Typical 
values  for  the  present  analysis  are 


kl  = O.Ol  and  k2  = 0.002 


(36) 


Finally,  k and  e are  computed  for  the  boundary-layer  region  from 


dv 

dy 


0.3  Uy./ 


(37) 


f = c - i- 
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5.1.2  Exit  Boundary  Condition 

With  a diffuser  incorporating  a constant  area  section  at  the  exit 
(tall  pipe),  the  tail  pipe  exit  boundary  condition  can  be  specified  by 
the  parallel  flow  condition,  i.e.,  3/3x  = 0, 

5.1.3  Symmetry  Condition 

The  symmetry  condition  is  imposed  to  reduce  the  computational 
domain.  Along  the  line  of  symmetry  (r  = 0),  the  following  conditions 
are  satisfied: 

W = 0 , 

= li  = 0 . 

dr  ^3 

dk  = df  ^ 

dr  dr 

V = 0 


5.1.4  Wall  Boundary  Condition 


Since  both  the  core  region  and  the  sublayer  region  are  solved 
numerically  with  the  present  formulation,  the  wall  boundary  condition  is 
applied  at  y = 0 as 


'V 


constant  , 


dx  dy 

= k = f = 0 , 
u = V = 0 


5.2  COORDINATE  TRANSFORMATION  WITH  A VARIABLE 
SUBLAYER  STRETCHING 


<39) 


In  the  present  analysis,  the  governing  equations  along  with  the 
boundary  conditions  are  solved  in  a transformed  rectangular  domain. 

This  is  achieved  by  transforming  the  diffuser  geometry  into  a rectangu- 
lar domain  followed  by  a sublayer  stretching  so  that  adequate  resolution 
can  be  obtained  in  the  wall  layer  region,  see  Fig.  3.  The  sublayer 
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Stretching  is  designed  to  adjust  automatically  for  the  sublayer  velocity 
gradient.  In  addition,  the  relative  number  of  grid  points  in. the  sub- 
layer and  the  core  region  can  also  be  varied  to  compensate  for  the  inlet 
boundary- layer  thickness. 

The  complete  coordinate  transformation  functions  in  the  core  and 
the  sublayer  regions  are  listed  below  along  with  the  transformation 
factors. 


Core  Region  Transformation  Functions 


(dll\ 

\dx/  fdr\  S (x)  \dr  / 

\d7) 


19 


AEDC-TR-77-78 


Sublayer  Region  Transformation  Function: 


r = S(x)  {1  I - tan/3 

;Tm7) 

I - -(ff)  • 

" « ■ S<*)  ■ (I;)  ■ 

-(§-:)  • {*"<^>  jfe  - • '} 


d^f  -s 
dx^ 


(41) 


j ^ * 2fi 

■ a • S (x)  ■ (li)  ®’"  [/5  ^"max- 


_SJx)  /^\ 
S(*)  \drJ 


where  a,  g,  c,  g and  F are  coefficients  of  transformation  (Ref.  1)  - 

These  coefficients  are  determined  through  a matching  procedure  at  location 

equal  to  unity  in  the  present  analysis).  By  selecting 

the  proper  value  of  r^^^  (which  must  be  greater  than  r^),  one  can  adjust 

the  relative  number  of  grid  points  in  the  sublayer  and  the  core  region. 

For  example,  with  r^  = 1 and  r^^^  = 2,  there  will  be  50  percent  of  the 

grid  points  in  both  the  sublayer  (r  < r < r ) and  the  core  (0  < r < 

« o — — max  — — 

r ) regions. 


6.0  NUMERICAL  METHODS 

6.1  A GENERAL  FINITE  DIFFERENCE  FORMULATION  WITH  DECAY 
FUNCTIONS 


The  governing  equations  in  the  transformed  plane  (x,  r)  can  be 
written  in  a standard  form  as 


(42) 
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where  <|>  represents  the  flow  variables  k,  and  e.  The  coefficients 

a^ , a^,  , b2,  and  d are  given  in  Appendix  A.  The  general  finite  dif- 

ference formulation  of  Eq.  (42)  is 


1 


(43) 


where  the  decay  functions  and  are  evaluated  locally  as 

C - 1.0  - 0.0625  (11)^  • |R|  ^2 


where  C 


2 

ml 


n\>2 


C.  • when  II 

3 


and  G 


G.  ^ when  R = 


_ 2 


8r 


(44) 


The  use  of  the  decay  function  assures  the  numerical  stability  as  well  as 
the  accuracy  throughout  the  flow  field  for  a wide  range  of  Reynolds 
numbers  (Refs.  8 and  9). 


6.2  GENERAL  PROCEDURE  WITH  A STANDARD  GAUSS-SEIDEL 
POINT  ITERATION  METHOD 


The  finite  difference  equation  [Eq.  (43)J  can  be  solved  In  several 
ways.  For  simplicity,  a standard  Gauss-Seldel  iteration  method  is  used 
to  obtain  the  solution.  The  method  updates  the  field  values  with  the 
most  recent  available  neighboring  point  data  in  the  following  way: 


n+l 

^1.) 


+ 


. n+l 


n+l  n 

d),  . ,+<!.. 

' J.j  - 1 i.j 


(45) 


cu 
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where  the  superscript  n denotes  the  values  at  nth  iteration,  (n+1) 
denotes  the  updated  value.  The  coefficients,  and  CU  are 


C 


2 


C 


3 


CU 


(46) 


The  general  flow  chart  which  describes  the  sequence  of  the  solution  pro- 
cedure is  presented  in  Fig.  4. 


7.0  FACTORS  AFFECTING  THE  RATE  OF  CONVERGENCE 
OF  THE  POINT  ITERATION  F1ETHOD 

The  advantage  of  using  the  Gauss-Seidel  point  iteration  method  to 
solve  coupled  algebraic  equations  is  well  known.  The  method  can  be  easily 
programmed  and  does  not  require  excessive  computer  storage.  However, 
there  are  two  Important  factors  which  affect  the  rate  of  convergence, 
namely,  (1)  the  total  number  of  grid  points  and  (2)  the  initial  flow- 
field  guess.  As  the  total  number  of  grid  points  is  increased  to  obtain 
better  resolution,  the  rate  of  convergence  decreases.  This  is  because 
the  information  at  the  boundary  points  must  travel  a longer  distance  to 
reach  the  whole  flow  field  in  the  elliptic  type  of  flow  calculations. 

Since  the  point  iteration  method  can  update  only  one  point  value  at  a 
time,  the  larger  the  number  of  grid  points  the  slower  the  convergence 
rate.  On  the  other  hand,  the  rate  of  convergence  can  be  greatly  im- 
proved with  a good  initial  flow-field  guess. 
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7.1  INITIAL  GUESS  FOR  DIFFUSER  FLOWS 


It  is  well  known  that  a turbulent  bound ary- layer  profile  can  be 
represented  by  the  superposition  of  a logarithmic  profile  and  a wake 
profile  (e.g..  Ref.  10).  The  representation  is  commonly  known  as  Coles' 
law  of  the  wake.  Coles'  concept  is  generalized  in  the  present  approach, 
with  the  logarithmic  part  of  the  profile  replaced  by  a function  related 
to  the  fully  developed  pipe  (or  channel)  flow  profile  at  the  diffuser 
inlet  and  with  the  wake  portion  represented  by  a cosine  function. 


The  initial  velocity  field  is  obtained  as  a linear  combination  of 
the  inlet  profile  and  the  wake  profile,  Fig.  5,  i.e., 


+ “b 


where  u^  is  related  to  the  inlet  velocity  profile,  u^,  as 
and  the  wake  component  of  the  profile,  is  defined  as 


(47) 


(48) 


(so)  “b.c 


I + cos 


(l-o-  ') 


(49) 


The  centerline  values  u^  ^ and  tij^  ^ are  the  only  unknowns  in  Eqs. 
(48)  and  (49).  From  Eq.  (47), 


a. c g b , c 


(50) 


With  Eq.  (50),  Eq.  (48)  becomes 


U — 11  , 

u = _£ . u 

" “l.c  ’ 


(51) 


Therefore,  the  final  expression  for  the  velocity  profile  can  be 
written  as 


^ - u.  ) 

u = . u,  + u 


I.C 


I “b.c 


1 + CDS  fJL.  . r\ 
\SO  / 


(52) 
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The  unknown,  ^ in  Eq.  (52)  can  be  determined  from  the  conservation  of 
mass.  By  integrating  Eq.  (52)  across  the  diffuser,  one  obtains 


u • rdr 


The  right-hand  side  can  be  evaluated  as 


SO(x) 


rdr 


SO(i) 


b , c 


[‘  * ■ ')1 


SO(x) 


rdr  + 


(53) 


rdr 


(5^) 


In  terms  of  the  stream  function,  4'  = (/urdr) , Eq.  (54)  becomes 

= (^^)-  (S0)‘  ■ {^)  d(-D* 

= AVj  (S0)2  ^ } 


(55) 


where 


1 

/U,  rdr 


Solving  for  one  obtains 


A’P 


4 


(56) 


^ » C 


When  the  wake  component  vanishes,  i.e.,  ^ = 0,  Eq.  (56)  becomes 


] - (SO)2  ^ = 0 


(57) 
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or 


(sor 


(58) 


In  general,  the  wake  component  (u,  ) is  nonzero  and  the  centerline 

D ) C 

velocity  (u^)  is  greater  than  that  given  by  Eq.  (58).  In  order  to 
increase  value  from  Eq.  (58),  an  effective  area  concept  is  introduced 
as 


U 

c 


(59) 


The  simplest  form  of  SO  ,,  is 

eff 


= [SO(k)  - 1)  - A -t-  1.0  (60)  * 

where  S0(x)  ,,  is  shown  in  Fig,  6 
ett 


The  unknown  constant  A (0  ^ A _<  1.0)  depends  upon  the  diffuser 
divergence  angle.  In  general,  the  value  of  A decreases  from  unity  as 
the  divergence  angle  increases.  No  formula  has  been  derived  to  deter- 
mine the  exact  value  of  A.  Therefore,  its  specification  requires  some 
experience  about  the  diffuser  flow  solution.  For  an  inexperienced  user, 
in  general,  A = 0,5  is  a good  Initial  guess. 


The  complete  initial  flow-field  distribution  of  the  velocity,  the 
vortlcity,  and  the  stream  function  (u,  0,  T)  can  be  derived  from  Eqs. 
(56)  and  (47)  as 


" (so)"  “1‘” 

a (jL)  . , 


-I- 


I.c 


'P,(r)  • 


u.  ) (S0)2 

+ 

Uf 

l,c 


sin 


(61) 

(62) 


(63) 
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(u  — u,  ) <1, 

£ ki£_  + 01* 

“l.c  ■ 2 


['  ■ (so  ■ ')] 


for  any  axial  location,  x. 


(64) 


7.2  MULTIPLE  GRID  GAUSS^EIDEL  ITERATION  METHOD 


The  Gauss-Seldel  iteration  method  is  a point  iteration  method  which 
updates  one  point  at  a time  with  the  most  recently  available  values  at 
the  neighboring  points.  It  can  be  easily  seen  that  the  rate  at  which 
the  information  from  a boundary  is  propagated  to  the  rest  of  the  flow 
field  depends  upon  the  total  number  of  the  grid  points  in  the  flow 
field.  When  the  flow-field  pattern  Is  fairly  complicated  and  it  is 
necessary  to  use  a large  niimber  of  grid  points  to  improve  the  accuracy  of 
the  final  solution,  the  rate  of  convergence  becomes  slow  with  the  point 
Gauss-Seidel  Iteration  method.  A simple  and  efficient  multiple  grid 
method  has  been  developed  to  speed  the  rate  of  convergence. 

A series  of  grid  systems  can  be  constructed  for  a rectangular 
domain,  Fig.  7,  with  grid  spacing  ratio  equal  to  (1/2)^, 

The  general  multiple  grid  Causs-Seldel  iteration  scheme  can  be 
written  as 

Do  100  J - JL  . .IR  . L 
Do  100  1 = IL  , IR  , L 


n+l 


n ii+1  n n+l 

4>,  . ^ {C.6.  , . + C.A.  , . + C,4>.  . , ' C,  <iS.  . , 

WJ  1 itL,;  2 L,j  3 ^i,j+L  4 ^i,j— L 

+ D")/(Cj  + Cj  + C3  + C^) 


(65) 


100  Continue 
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With  L equal  to  unity,  one  recovers  the  conventional  single  grid  system 
Gauss-Seidel  method.  For  a dual  system,  L is  set  equal  to  two  for  the 
first  N iterations  to  obtain  a solution  for  the  coarse  grid  system. 
Before  continuing  to  a finer  grid  system  (with  L = 1),  it  is  necessary 
to  obtain  the  midpoint  value  with  a linear  Interpolation  scheme.  The 
flow  chart  for  the  process  is  shown  in  Fig,  8.  Since  the  coarse  grid 
system  has  only  one-fourth  of  the  total  grid  points  of  the  finer  system, 
the  rate  of  convergence  is  a factor  of  four  faster.  Another  advantage 
is  that  the  fast  convergent  solution  obtained  in  the  coarse  grid  system 
is  used  to  obtain  the  initial  guess  for  the  finer  grid  system  through  an 
interpolation  routine.  In  the  present  solution  procedure,  a dual  grid 
system,  i.e.,  from  a 1 7 by  31  coarse  system  to  a 33  by  61  fine  grid 
system  is  adequate.  The  operation  is  done  automatically  in  the  computer 
program. 


8.0  RESULTS  AND  DISCUSSION 

The  numerical  formulation  and  the  turbulence  model  presented  in 
the  previous  sections  are  applied  to  the  computation  of  several  conical 
diffuser  flows. 

Although  much  experimental  data  on  the  performance  of  subsonic 
diffusers  are  available,  few  of  these  data  are  suitable  for  verification 
of  the  numerical  analysis.  The  major  emphasis  in  most  experimental 
investigations  has  been  placed  on  overall  performance,  i.e.,  pressure 
recovery,  and  only  rarely  were  the  mean  velocity  profiles  at  the  inlet 
section  and  in  the  divergent  section  measured.  Even  more  rare  are 
experiments  in  which  the  turbulence  properties  in  the  flow  field  were 
measured.  Recently,  however,  experimental  data  on  both  the  mean  velocity 
and  turbulence  properties  in  an  8-deg  conical  diffuser  were  obtained  by 
Okwuobi  and  Azad  (Ref.  4);  the  inlet  flow  was  a fully  developed  pipe 
flow.  These  data  provide  a good  test  case  for  comparison  with  the 
numerical  analysis  because  the  fully  developed  entry  conditions  are  well 
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established  and  can  be  predicted  with  confidence.  In  addition  to  the 
correlations  of  the  Okwuobi  and  Azad  experiment,  comparison  of  the 
numerical  results  with  other  experiments  (in  which  the  inlet  conditions 
are  less  well-defined)  are  described  in  this  section. 

8.1  NUMERICAL  SOLUTIONS  OF  8-DEG  CONICAL  DIFFUSER  WITH  A 
FULLY  DEVELOPED  INLET  PROFILE 


The  geometry  of  the  8-deg  conical  diffuser  is  shown  In  Fig.  9. 

The  area  ratio  is  4:1.  The  computed  velocity  field  at  several  axial 
stations  is  also  shown  in  Fig.  9.  It  can  be  seen  that  the  velocity 
profile  changes  from  that  of  a fully  developed  inlet  pipe  flow  into  a 
wake  profile  in  the  diverging  section.  The  law  of  the  wall  region  also 
diminishes  and  the  wake  portion  of  the  profile  becomes  dominant  at  the 
exit  section  of  the  diffuser.  The  comparison  between  the  computed  and 
the  measured  velocity  profiles  at  two  stations  is  shown  in  Fig.  10.  The 
agreement  between  the  experiment  and  theory  is  excellent.  In  order  to 
check  the  sensitivity  of  the  computed  flow  field  to  the  present  turbu- 
lence model,  a small  perturbation  on  the  c formula  from 


to 


= c - 1.4  Ac 

/I 


(66) 


was  used.  The  modification  has  no  effect  on  the  fully  developed  inlet 
velocity  profile  because  the  flow  is  parallel  so  that  the  radial  veloc- 
ity component,  v,  vanishes.  The  velocity  profile  at  the  downstream 
station  (x/Dj  = 5.95)  is  only  slightly  changed' by  the  c^  perturbation  as 
shown  in  Fig.  10. 


The  velocity  profile  near  the  wall  is  shown  in  Fig.  11.  The  line- 
arity of  the  velocity  profile  in  the  sublayer  region  Is  clearly  evident 
and  the  agreement  with  the  experimental  data  is  excellent. 
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The  decay  of  the  centerline  velocity  is  shown  in  Fig.  12  for  two 
models.  The  = c^^  model  gives  a faster  velocity  decay  as  expected. 

On  the  other  hand,  results  obtained  from  the  new  c^  expression  agree 
very  well  with  the  experimental  data. 

The  total  shear  stress  distribution  is  given  in  Fig.  13.  The  com- 
puted results  agree  well  with  the  measured  total  shear  stress  at  x/D^.  = 
1.19  in  the  diverging  section.  At  the  downstream  station  (x/D^.  = 

5.95),  although  the  shape  of  the  total  shear  stress  profile  is  qualita- 
tively correct,  the  maximum  predicted  value  is  about  15  percent  higher 
than  the  measured  value.  The  measured  turbulent  kinetic  energy  dis- 
tribution, Fig,  14,  is  also  lower  than  the  prediction  at  the  downstream 
station,  especially  near  the  wall.  In  this  region,  the  turbulent  fluctua- 
tion velocities  are  of  the  same  order  of  magnitude  as  the  local  velocity 
(u).  In  this  situation,  the  hot-wire  probe  measurement  is  questionable. 
Usually,  the  x-probe  yields  a lower  measurement  of  the  turbulent  intensity 
because  of  the  higher  instantaneous  flow  angularity  and  negative  velocity, 
hence  an  underestimated  shear  stress  and  TKE.  The  agreement  between 
calculated  and  measured  TKE  near  the  diffuser  inlet  is  quite  good.  The 
calculated  turbulent  kinetic  energy  distribution  near  the  wall,  shown  in 
Fig.  15,  is  lower  at  the  downstream  station  than  at  the  inlet.  This  is 
caused  by  the  lower  TKE  production  associated  with  the  downstream  veloc- 
ity profile.  The  ratio  between  the  shear  stress  u*v*  and  the  TKE  1s 
shown  in  Fig.  16,  At  the  inlet,  the  ratio  quickly  rises  to  a maximum 
value  of  0,3  and  drops  gradually  to  zero  at  the  centerline.  Experimental 
data  for  fully  developed  pipe  flow  by  Laufer  (Ref.  11)  are  also  Included 
in  Fig,  16  for  comparison.  At  the  downstream  station,  the  profile  is 
skewed  toward  the  centerline  and  the  maximum  value  of  the  ratio  also 
increases  slightly. 

In  Fig.  17,  the  distribution  of  the  eddy  viscosity  coefficient, 
c , is  shown  for  several  axial  stations.  At  the  inlet,  the  c 

M U 

distribution  is  identical  to  that  of  the  previous  model  since  the  Ac 
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term  is  zero.  The  distribution  drops  immediately  in  the  diverging 
section.  The  level  of  then  increases  slowly  with  increasing  distance 
in  the  diverging  section.  It  is  important  to  note  that  the  average 
value  of  c^  in  the  diverging  section  is  always  lower  than  the  inlet  c^ 
level. 

The  wall  pressure  coefficient,  C^,  is  presented  in  Fig,  18.  The 
computed  pressure  coefficients  agree  well  with  measured  values  except 
that  they  are  somewhat  overpredicted  in  the  downstream  stations.  Part 
of  this  difference  could  be  attributed  to  the  effect  of  the  computa- 
tional grid  size.  Another  factor  is  that  the  numerical  computation  is 
carried  out  for  a diffuser  with  a tail  pipe  while  the  experimental 
configuration  did  not  have  a tail  pipe.  It  has  been  shown,  Refs.  12 
and  13,  that  the  tall  pipe  Increases  the  pressure  recovery  somewhat. 
Therefore,  the  higher  predicted  pressure  recovery  is  consistent  with 
experimental  evidence. 

8.2  NUMERICAL  SOLUTION  OF  8-DEG  CONICAL  DIFFUSER  WITH  A 
BOUNDARY-LAYER  INLET  PROFILE 

Computations  of  an  8-deg  conical  diffuser  flow  field  were  made  with 
two  inlet  boundary-layer  thicknesses.  The  velocity  distribution  for  the 
thin  boundary-layer  case,  6^  = 0.05  D^.,  is  shovm  in  Fig.  19.  The  dif- 
ference between  the  velocity  profile  with  a boundary  layer  and  the 
fully  developed  one  is  the  existence  of  a potential  core.  It  can  be 
seen  from  Fig.  19  that  the  potential  core  remains  throughout  the  dif- 
fuser section.  The  wall  pressure  coefficient,  C^,  is  shown  in  Fig.  20. 
Comparison  with  the  experimental  pressure  recovery  is  good. 

The  velocity  distributions  for  the  thicker  boundary-layer  case, 

(6j.  = 0.2  Dj)  are  shown  in  Fig,  21.  The  potential  core  region  is 
smaller  than  in  the  thin  boundary- layer  case,  otherwise,  the  development 
of  boundary-layer  profiles  is  about  the  same.  The  total  shear  stress 
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distribution  is  shown  in  Fig.  22.  Upstream  of  the  inlet,  the  total 

shear  stress  is  highest  and  constant  near  the  wall,  l.e..  t = t . 

tot  w 

In  the  diffuser  section,  the  maximum  total  shear  stress  location  moves 
away  from  the  wall  and  the  maximum  value  increases.  The  shear  stress 
along  the  wall  decreases  in  the  diverging  section.  The  turbulent 
kinetic  energy  distribution  is  shown  in  Fig.  23.  The  free-stream 
turbulence  level  was  set  equal  to  one  percent  which  remains  approxi- 
mately constant  in  the  potential  core  region.  Just  as  was  observed  for 
the  shear  stress,  the  location  of  the  maximum  TK£  moves  away  from  the 
wall  in  the  diverging  section.  The  eddy  viscosity  distribution,  nor- 
malized by  the  local  centerline  velocity  and  the  local  boundary- layer 
thickness,  is  shown  in  Fig.  24.  Notice  that,  as  x Increases,  the  maximum 
eddy  viscosity  increases  greatly  from  the  boundary-layer  value  of  order 
0.001  near  the  entrance  to  almost  the  turbulent  mixing  layer  value  of 
order  0.01  at  the  downstream  station.  The  length  scale  derived  from 
the  turbulent  kinetic  energy  dissipation  is  shown  in  Fig.  25.  The 
maximum  normalized  length  scale  in  the  boundary  layer  increases  from 
0,65  to  about  1.4  between  the  inlet  station  and  the  downstream  station. 
Finally,  the  wall  pressure  coefficient,  C^,  is  presented  in  Fig.  26.  The 
exit  value  of  0.77  is  slightly  lower  than  the  previous  case  with  a 

thinner  inlet  boundary  layer  (C  = 0.8).  This  is  consistent  with 

P 

experimental  evidence  that  decreases  with  increasing  inlet  boundary- 
layer  thickness  and  that  increases  again  when  the  inlet  condition  is 
fully-developed  (Refs.  13  and  14).  The  effect  of  the  inlet  boundary- 
layer  thickness  on  the  pressure  recovery  coefficient  for  an  8-deg 
conical  diffuser  is  shown  in  Fig.  27.  Although  the  trend  is  well  pre- 
dicted, the  predicted  level  Is  higher  than  the  data  of  Ref.  13. 

8.3  NUMERICAL  SOLUTION  OF  16-DEG  CONICAL  DIFFUSER  WITH  A 
BOUNDARY-LAYER  INLET  PROFILE 

The  computed  velocity  distribution  for  a 16-deg  conical  diffuser 
with  a thin  inlet  boundary  layer  (5^  = 0.05  D^.)  is  shown  in  Fig.  28. 
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The  velocity  at  the  exit  station  (x/D^  = 3.58)  is  negative  near  the  wall 

indicating  that  the  flow  is  separated.  This  is  in  agreement  with  the 

experimental  evidence  observed  by  McDonald  and  Fox»  Ref.  15,  and  Senoo 

and  Mishi,  Ref.  (12),  who  found  the  16-deg  conical  diffuser  with  an  area 

ratio  of  four  lies  within  the  separated  flow  regime  of  the  diffuser 

performance  map.  The  pressure  coefficient  C is  shown  in  Fig.  29. 

P 

The  agreement  with  the  experimental  data  (Ref.  12)  is  qualitatively 
good.  The  experimentally  observed  separation  point  detected  by  tufts 
lies  somewhat  ahead  of  the  predicted  point  of  separation.  This  could  be 
caused  by  the  coarse  grid  size  (Ax)  used  in  the  computation.  However, 
the  use  of  tufts  only  provides  a crude  estimate  of  the  point  of  separation. 
Thus,  the  agreement  is  probably  as  good  as  one  could  expect. 

8A  NUMERICAL  SOLUTION  OF  13-DEG  CONICAL  DIFFUSER  WITH  A 
BOUNDARY-LAYER  INLET  CONDITION 


The  measured  and  calculated  velocity  distribution  for  a 13-deg 
conical  diffuser  is  shown  in  Fig.  30.  The  agreement  with  the  data. 

Ref.  16,  is  qualitatively  good  at  x/D^  - 3.2.  The  experimental  data, 
obtained  with  an  x-probe  hot  wire,  shows  a larger  tail  region  near  the 
edge  of  the  boundary  layer.  Near  the  wall,  the  hot-wire  measurement  is 
slightly  higher  than  the  calculated  velocity  values  because  of  high 
local  relative  turbulent  intensity.  At  a far  downstream  station,  x/D^ 

= 5.9,  the  flow  separation  is  clearly  indicated  in  the  calculated  velocity 
profile.  Under  the  separated  flow  condition,  the  conventional  hot- 
wire probe  is  not  capable  of  distinguishing  the  positive  and  negative 
velocities.  As  a result,  the  hot-wire  measurement  is  always  positive  in 
the  separated  flow  region.  The  pressure  coefficient  is  shown  in  Fig. 

31 . The  agreement  between  the  data  and  the  computed  value  is  good.  The 
predicted  value  of  is  slightly  higher  than  the  data  at  the  exit.  The 
point  of  separation  is  also  Included  in  Fig.  31.  The  present  predicted 
point  of  separation  lies  among  the  separation  points  predicted  with 
various  assumptions  in  the  boundary-layer  theory  (Ref.  16).  No  experi- 
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mental  evidence  of  the  point  of  separation  was  obtained.  Furthermore, 
no  comparisons  were  made  between  the  computed  turbulence  flow  field  and 
the  data  because  of  the  apparent  misalignment  problem  of  the  hot-wire  x- 
probe  used  in  the  experiment  as  discussed  in  Boldman's  report  (Ref.  16). 

9.0  CONCLUDING  REMARKS 

The  important  results  of  this  analytical  and  numerical  study  of 
turbulent  flows  in  conical  subsonic  diffusers  are: 

1 . A c^  formula  for  axisymmetrical  flows  has  been  developed  for 
the  low  Reynolds  number,  two-equation  k-e  turbulence  model 
which  unifies  the  computation  of  2D  planar  and  axlssrrametri- 
cal  flows, 

2.  An  initial  guess  procedure  and  a multiple  grid  Gauss-Seldel 
iteration  method  have  been  developed  to  increase  the  rate  of 
convergence. 

3.  Numerical  computations  of  several  separated  and  nonseparated 
conical  diffuser  flow  fields,  with  nonuniform  inlet  condition 
varying  from  a fully  developed  pipe  flow  profile  to  a thin 
boundary-layer  profile,  are  in  good  agreement  with  available 
experimental  data. 
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Figure  1.  Centerline  velocity  distribution  in  8-deg  conical  diffuser. 
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Figure  2.  Distribution  of  the  eddy  viscosity  coefficient,  for  a round  jet. 
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Figure  4.  Flow  chart  of  the  solution  procedures. 
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Figure  6.  Effective  area  concept  for  initial  centerline 
velocity  distribution. 
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Figure  8.  A dual  grid  Gauss-Saidel  itaration  schame. 
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Figure  11.  Velocity  distributions  near  the  wall  in  8-deg  conical  diffuser. 
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Figure  13.  Total  shear  stress  distributions  in  8-deg  conical  diffuser. 
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Figure  14.  Turbulent  kinetic  energy  distributions  in  8-deg 
conical  diffuser. 
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Figure  15.  Turbulent  kinetic  energy  distributions  near  the  wall 
in  8-day  conical  diffuser. 


AEDC-TR-77-78 


Turbulent  Shear  Stress 
TKE 


Figure  16.  Turbulent  shear  stress-TKE  ratio  in  8-deg 
conical  diffuser. 
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Figure  19.  Velocity  distribution  in  8-deg  conical  diffuser  with  a 
boundary-layer  inlet  profile,  5|  = 0.05  D| . 
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Figure  20.  Pressure  coefficient  in  8-deg  conical  diffuser 

with  a boundary-layer  inlet  profile,  5|  - 0.05  D| . 
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Figure  21 . Velocity  distributions  in  8-deg  conical  diffuser  with 
a boundary-layer  inlet  profile,  5|  = 0.2  D| . 
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Figure  22.  Total  shear  stress  distributions  in  8-deg  conical  diffuser 
with  a boundary-layer  inlet  profile,  5|  = 0.2  D| . 
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Figure  24.  Eddy  viscotity  distributions  in  8<teg  conical  diffuser  with 
a boundary-layer  inlet  condition,  5|  - 0.2  D|. 
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Figure  25.  Length  scale  distributions  in  8-deg  conical 

diffuser  with  a boundary-layer  inlet  condition, 
5i  = 0.2  D,. 
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Figure  26.  Pressure  coefficient  in  8-deg  conical  diffuser  with  a 
boundary-layer  inlet  profile,  5|  = 0.2  D|. 
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Fi^re  27.  Effsct  of  inlet  boundary -layer  thickness  on  the  pressure 
recovery  coefficient  for  8-deg  conical  diffuser. 
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Figure  28.  Velocity  distributions  in  16-deg  conical  diffuser  with 
a boundary-layer  inlet  profile. 
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Figure  29.  Pressure  coefficient  in  l6-deg  conical  diffuser 
with  a boundary-layBr  inlet  condition. 
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Figure  30.  Velocity  distributions  in  13-deg  conical  diffuser 
with  a boundary-layer  inlet  condition. 
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Figure  31.  Pressure  coefficient  for  a 13-deg  conical  diffuser  with 
a boundary-iayer  iniet  profile. 
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APPENDIX  A 

GOVERNING  EQUATIONS  IN  TRANSFORMED  COORDINATES 


The  standard  form  of  the  derived  equations.  Ref.  1,  is 

1“' 5?  " i'"' H air  ° 


<A-1) 


where  (ji  represents  the  flow  variables,  i.e.,  fJ,  S',  k,  and  e.  The  cor- 
responding coefficients,  , a^,  b^,  b^,  and  d are  given  below  for  each 
of  variable,  n,  ¥,  k,  and  c. 

(^-Equation 

ai  = I 


1 (u 


^ [w  * (ft)l^]](f^)-  [(S) " (&0]| 


, ^ /dr  ^ dv  1 

(v--Vj)  ^ ,Jx^  Ldi  \5r/  d»  \dx  / df  J 


(A-2) 


* = (l^)(li)[fF(if)-  (ii)-(i^)  i?]} 

2 

if d^f  _ dfF\  rdu  /dr  \ + dv  . /dr  \ dv 

P2  lbx2  Ur^  ^ id7>> 

’^  (&)[!?  (if)  -(fi)-(f^)|?]} 

(If) [i? (if)- [ft ^ (£)  fr]]}( 


+ 2 
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V-Equation 


a,  - 1 


■[(E)  - (If)] 


b,  - 0 


b = ~ - ill  + ^ f iL\ 

- ■ ^ 1^- (E)  * 


(A-3) 


k-Equation 


^ = 1 


(i'  + 


•.-[(E)“^(EX] 

* (E)  1^]} 

1.7*77  |[»-(E)^-^--p](E) 

^["-(§^^(E)l*)](E)-(-^.)[(S)*(0)]| 
^ (E)  * K(E  ^ (E)  Ey  * [E  (E)]’  * * (01 

+ 4-  ^ + - 1 /ifi^  + ^ I 1, 

L^f  dx  \dx  / 3rJ  t (i/ + i/j)  (v  + Vj)  \ y2  j 


(A-4^) 


c-Equation 


aj  = 1 


- [(E)’  * (E)’] 

7^{"-[A(^)^(E)E(^)]} 
w|['-(E)E(^)-f‘^-v<^p](E) 
^[“-(E(^)*(E)E(^))](E) 

- ‘-•'.'V[(|:o^(0)]| 


'>2  ■ 7 

2 iv  ■>!■ 


(A-5) 
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d = 2 / + c I2  r^—  +■  ( 

\()\/  dxdf  ' (i-' + ■'crl  I LV^ji  \ci\/  d’/ 

*fe(ff)r  * ® (f)’]^  [i?(i)*  fi* 

“C  ^ I 

2 k (t/  + i/|  . (j^) 

In  the  vorticity-stream  function  formulation,  the  pressure  field  is 
decoupled  from  the  momentum  equations.  Therefore,  the  pressure  field 
can  be  calculated  separately  from  the  pressure  equation.  In  the  trans- 
formed plane,  the  pressure  equation  can  be  vnc itten  as 


i * ■ m * (ffn 


& - s 


(^)  (If  n 

+ 2 ^ - /v^  + u^N  / 

\dx  / dxdr  \ 2 / j 


(A-7) 


5u^  + v|^+ 

dr  V dr  / r dx  Vdx  / df 

The  corresponding  coefficients  for  the  pressure  equation  in  the  standard 
form  are: 

“1  = ^ 

b,  = 0 

^=--[(0)^(f?)*Hlf)]  (A-S) 

• ^ “ (If ) Ifc-  - * j|  (^)-  # (^“)  ■ [(lf)“  ^ (If)’] 

■ * i ("V^)[(fr() " (0) " ^ (if)] * ^ (If)  & (^)j 
-»f(^)-S“^'i-(lf)f 
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^1’  ^2 
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Ci,C2 


C • C ■ c ^ 
p’  yo’  y1 

Ac 
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Gi.Gj 


NOMENCLATURE 

Parameter 

Constant 

Coefficients 

Constant 

Coefficients 

Transformation  coefficient 
Coefficients 

Pressure  Coefficient,  — — - 

0.5  pu^ 

Eddy  viscosity  coefficients 

Eddy  viscosity  coefficient  correction  for  axlsyrmtietrlc 
flow 

Local  diffuser  diameter 
Source  term 

Transformation  parameter 

Rodi's  c^  correction  for  axlsymmetrlc  jet 

Decay  Functions 
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g 

K 

k 

kl,  k2 


I 


n 


P 

R 

Re 


R^.Rj 


r 

max 


S(x),  SO(x) 


u 


u,v 


Au 


Transformation  parameter 

Constant  In  Schllchtlng ' s round  jet  theory 

Turbulent  kinetic  energy 

Proportionality  constants  for  k and  In  potential  core 

Length  scale 

Iteration  number 

Pressure 

Parameter 

DjUi 

Reynolds  number,  — — 

Grid  Reynolds  number 

Radial  coordinate  in  physical  and  transformed  plane, 
respectively 

Diffuser  wall  coordinate  in  transformed  plane 

Diffuser  wall  coordinate  in  physical  plane 

Average  velocity 

Velocity  components 

Velocity  excess  in  Rodi's  c^  model 
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U 


.W 


u 


V 


b 

* 


X 


y 

YG 


a, 6 


6 


e 


26 


V 


c 

e 


T 

♦ 


Turbulent  velocity  components 
Wake  component  of  the  velocity  profile 
H Friction  Velocity 

Axial  coordinate 

Lateral  distance  measured  from  wall 

The  width  of  a round  jet  in  Rodl's  c model 

Transformation  parameters 

Index,  zero  for  planar  configuration,  1 for  axlsymmetric 
configurations,  or  incremental  or  boundary- layer  thickness 

Turbulent  kinetic  energy  dissipation 

Total  diffuser  divergence  angle 

Molecular  viscosity 

Eddy  viscosity 

Constant 

Shear  stress 

Dependent  variable 

Stream  function 
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Q 

n 

p 

X 

SUBSCRIPTS 

c 

e 

I 

max 

o 

T,  tot 
w 

SUPERSCRIPTS 


Vortlclty 

Similarity  parameter  for  round  jet 

Density 

Constant 

Centerline 

Edge  of  the  boundary  layer 

Inlet  section 

Maximum 

Sublayer^  core  region  solution  matching  location 

Total 

Wall 


+ 


Dimensionless  quantity  with  the  sublayer  scale 


